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Structural, electronic and magnetic properties were calculated for the optimized a-U/W(110) thin 
films (TFs) within the density functional theory. Our optimization for 1U/7W(110) shows that the 
U-W vertical interlayer spacing (ILS) is expanded by 14.0% compared to our calculated bulk W-W 
ILS. Spin and orbital magnetic moments (MMs) per U atom were found to be enhanced from zero 
for the bulk of a-U to 1.4 \ib and -0.4 \ib at the interface of the 1U/7W(110), respectively. Inversely, 
our result for 3U/7W(110) TFs shows that the surface U-U ILS is contracted by 15.7% compared to 
our obtained bulk U-U spacing. The enhanced spin and orbital MMs in the 1U/7W(110) were then 
found to be suppressed in 3U/7W(110) to their ignorable bulk values. The calculated density of 
states (DOS) corroborates the enhancement and suppression of the MMs and shows that the total 
DOS, in agreement with experiment, is dominated in the vicinity of Fermi level by the 5f U states. 
Proximity and mismatch effects of the nonmagnetic W(110) substrate were assessed and found to 
be important for this system. 

PACS numbers: 73.20.At, 68.35.bd, 75.70.Ak, 75.70.Rf, 71. 15. Mb, 71.27. +a 



I. INTRODUCTION 

Uranium, as the most predominant light actinide, 
has attracted considerable interests experimentallyir— 
and thcorcticallj«^£i2r— in both of its bulk£~— and 
surfaced— ^ states due to the dual nature^ of its 5f elec- 
trons. An experimentally detected sharp peak at 2.3 eV 
below the Fermi energy (Ep) for a thin foil of a-U re- 
vealed that the 5f electrons could be nearly localized at 
the surface of this system*. However, it was very soon 
realized that the peak was originated from the 6d valence 
band emission due to the exposure of uranium surfaces 
to O2 which conversely implied an itinerant behavior for 
the 5f-electrons at the U surface^. A dispersive behavior 
was also foundi for the U 5f states and attributed to the 
direct f-f interaction, which excellently verified the band- 
like feature of the 5f U states reported in Ref^. Fur- 
thermore, the first preformed scanning tunneling spec- 
troscopy experiment to observe 5f states^ for the evap- 
orated well ordered uranium films up to 80 A on the 
W(110) crystal at room temperature reconfirmed that 
the 5f electrons would not be considered as localized elec- 
trons. Therefor, according to the aforementioned exper- 
imental evidences^£ we do not treat a-U/W(110) thin 
films (TFs) as a strongly correlated system in this work, 
which effectively saves the time of calculations in our 
ab initio surface calculations. There are several exper- 
imental evidences^— which reveal that the structure of 
U ovcrlayers is hypothetical hep in the U/W(110) TFs. 
On the other hand, there are some contradicting exper- 
imental results^ which show that the structure of U 
overlayers is orthorhombic. The hcp-U TFs&£ and a- 
U TFsii in their hypothetical freestanding forms have 
been theoretically studied. The orthorhombic a-phase is 
more stable than the hep- and 7-phases^. Springell et 
al. within the full-potential (FP) linear muffin-tin or- 



bital (LMTO) method calculated spin magnetic moment 
(MM) to be ~0.1/is^ for the bulk of the 7-U compound 
and to be more enhanced to for the surface of 

the freestanding 7-U TFs. Stojic and coworkers employ- 
ing the FP-augmented plane waves (FP-APW) method 
have shown that the a-U compound is a nonmagnetic 
system in its bulk form, whereas it behaves as a mag- 
netic system with the imposed spin MM of 0.84//^ at its 
freestanding surfaces^. A typical substrate can affect 
the structure^ and magnetism^ of the overlayers that 
can give rise to a variety of lattice mismatches and as a 
result to various U-U interlayer (IL) spacings (ILSs). The 
U-U spacing can change the degree of 5f localizations and 
thereby electronic structures. By going beyond the hy- 
pothetical freestanding hep- and a-U surfaces, we then 
aim to systematically study the surface properties and 
electronic structures of the realistic a-U/W(110) TFs. 
On the contrary to the freestanding hep- and a-U TFs, 
our result shows that the 3U/7W(110) is not a magnetic 
system. 

Our work is organized as follows. We first describe 
computational details in Sec. [TI] Subsequently in Sec. IIIII 
we elaborate on the simulation procedure. This section is 
presented together with the structural properties of the 
W and a-U in their bulk and surface states to optimize 
the superlattice structure of the a-U/W(110) TFs. In 
Sec. IIVI changes of ILSs due to the relaxation in differ- 
ent slabs are calculated and discussed. We work out the 
magnetic properties by shedding light into the electronic 
structures of the a-U/W(110) TF s in Sec. El which is 
divided to 5 subsections. In Sec. IV Al we focus on the 
bulks of U and W as the fundamental materials of this 
surface study. In Sec. IV Bl we concentrate on the hypo- 
thetical freestanding a-U surfaces to somehow validate 
our subsequent predictions on the actual a-U/W(110) 
thin solid films (TSFs). In this section we in excellent 
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agreement with the reported results in Reb^i find out 
that the hypothetical freestanding a-U(OOl) can behave 
as a magnetic system. In Sec. IV CI we shall turn our at- 
tention to a-U/W(110) TSFs, as the main subject of this 
work. In Sec. lVCl we show that the a-U/W(110) is not a 
magnetic system, which is in contrast to the freestanding 
result. The source of such a different magnetic behavior 
is the subject of Sec. IV Dl where the hypothetical free- 
standing a-U(OOl) and the actual a-U/W(110) TSFs are 
compared. In this section we concentrate on the verti- 
cal IL U-U distances as the main reason for the different 
behavior. We then complete the work in Sec. IV El with 
a discussion on the effects of tungsten substrate on its 
uranium overlayers. In this section we show the effect of 
tungsten on the in-plane U-U distances and as a result on 
the magnetization of the a-U/W(110) system compared 
to that of the freestanding a-U(OOl) TFs. Finally, we 
summarize and conclude the work in Sec. I VII 



II. COMPUTATIONAL DETAILS 

The FP-APW + local orbitabi^ (FP-APW+lo) 
method is used for solving the set of density func- 
tional theory (DFT) equation s 17 ' 18 as embodied in the 
WIEN2k codeiS to solve the Kohn-Sham equations^. 
For the exchange-correlation potential no more corre- 
lations are included than what have been implemented 
in the Perdew and Wang local density approximation^ 
(PW-LDA) and the Perdew-Burke-Ernzerhof generalized 
gradient approximation^ (PBE-GGA), because of the 
reported band-like character in Refsj2~— for the 5f U 
electrons in U/W(110) system. The orbital polarization 
has been estimated using a non-self-consistent treatment, 
which may be reasonable due to the assumed itinerant 
character for the 5f states. The relativistic spin-orbit 
(SO) interactions were incorporated to the band struc- 
ture calculations in a second variational procedure. The 
Muffin-tin radii were adjusted to Rmt=2.4 bohr for the 
W and to Ra/t=2.5 bohr for the U atoms. The expan- 
sion of the wave functions and charge densities were cut 
off by the RjifTKmai = 7 and G max = 12 parameters, re- 
spectively. The full relaxations were performed with the 
criterion of 1 mRy/bohr on the exerted forces. A set of 
14x14x1 special k-points has been used for integrations 
over the Brillouin zone of the lxl surface cell. 



III. OPTIMIZATION OF SUPERLATTICE 
SURFACE STRUCTURE 

We calculated the lattice parameter of the bcc-W to be 
3.154 A within LDA and 3.176 A within GGA. These re- 
sults are in agreement with the pseudopotential results of 
the others^, i.e., 3.129 A within LDA and 3.175 A within 
GGA. The GGA results arc in more agreement with the 
experimental valued of 3.165 A than those of the LDA. 
We cut the bec unit cell of the tungsten by the (110) 
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FIG. 1: (Color online) Black and white circles represent tung- 
sten (W) and uranium (U) atoms, respectively, (a) Crossview 
of bec-tungsten along [110] direction as our clean W(110) sub- 
strate with a = c = y/2 a,bulk = 4.492 A and b = bbuik = 
3.176 A. (b) Orthorhombic structure of a-U with a = 2.842 A, 
b = 5.780 A, c = 4.985 A and the relaxed internal parameter 
y = 0.105. (c) Doubled superlattice structure of 3U/7W(110) 
TFs by depositing a-U after shifting its middle plane to the 
right-hand side by the (1/2 - 2y)b on W(110) taking lattice 
mismatch and enough amount of symmetric vacuum, 20 bohr 
= 10.58 A, into account. Left highlighted portion constitutes 
the conventional unit cell of our Q-U/7W slab. The primi- 
tive unit cell of the highlighted slab constitutes our superlat- 
tice structure. The super lattice structure crystallizes in the 
Cmmm space group and contains totally 13 atoms, 2x3=6 
U atoms and 2x3+1=7 W atoms. U(W) represents image of 
U(W) with respect to an imaginary central plane which hori- 
zontally crosses the slab. Zo (Z0.5) shows the midpoint of the 
substrate (vacuum) of the 3U/7W(110) system. 



plane and reconstructed another unit cell along the [110] 
crystallography direction. As shown in Fig.[T](a), the re- 
constructed unit cell as our clean W(110) substrate looks 
like the fee structure with a = c / b. Lattice parameters 
of the fcc-like W(110) structure were fixed to be a = c = 
V% &buik = 4.492 A and b = hbuik = 3.176 A, where our 
GGA lattice parameter for the bcc-W, 3.176 A, was used 
for &buik and bt, u ik- Seven monolayer (ML) W(110) films 
were found quite enough to model the substrate. We fully 
relaxed the orthorhombic a-U, as shown in Fig.Q](b), and 
calculated the lattice parameters including internal pa- 
rameter y to be (a, b, c, y) = (2.842 A, 5.780 A, 4.985 A, 
0.105) in agreement with the experimental values^ of 
(a, b, c, y) = (2.844 A, 5.869 A, 4.932 A, 0.102). Our 
calculated equilibrium volume of the primitive unit cell 
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for the a-U, V Ua = abc/4 = 20.472 A 3 , is in excellent 
agreement with the experimental valued, 82.3256/4 = 
20.581 A 3 , and can be compared with the values ob- 
tained from other theoretical methods^— as well. The 
a-U structure represented in Fig. [T] (b) can be trans- 
formed to the fcc-like structure, if its middle horizontal 
plane is shifted to the right-hand side along y-axis by 
the (1/2 - 2y)b. In this case, the Fig. Q] (b) looks like 
Fig. Q] (a), as both of them are now fcc-like. However, 
their lattice parameters are not identical to each other, 
as discussed above and shown in Figs. [TJ Therefore, the 
substrate, due to its different in-plane and vertical lattice 
parameters compared to those of the U overlayers, gives 
rise to lattice mismatches. By considering the in-plane 
lattice parameters, a = 4.492 A, b = 3.176 A, and verti- 
cal ILS, lattice parameter, c/2 = 4.492/2 A = 2.246 A, 
of the clean W(110) surface, 2x3 ML of the a-U were 
symmetrically deposited on top and bottom of the 7 ML 
W(110). The three U layers on bottom of the 7W(110), 
i.e. Ul, U2, U3, which can be considered as the images 
of those on top of the substrate, were also deposited to 
add inversion symmetry. In order to reach the highest 
possible symmetry, less traditional slab is set up. The 
vacuum is then added in the interior of the slab, as shown 
in Fig. [T] (c). In traditional slab the vacuum is usually 
added on top and below of the slab. Here, in contrast 
to the traditional slab, the vacuum is created in between 
of the outer U layer (U3) and its image U3, as shown in 
Fig. [T] (c). The Wi also represents the image of Wi with 
respect to an imaginary central plane which horizontally 
crosses the slab for i=l, 2, 3, 4. The vacuum thickness 
was determined to be 10.58 A by calculating the total 
energy, work function, and exerted forces on the surface 
atoms versus various vacuum thicknesses for the 7 ML 
W(110) clean substrate. The total energy, work function 
and forces were well converged at the determined vacuum 
thickness - they did not considerably fluctuate by increas- 
ing the vacuum thickness more than the value obtained. 
The jU/iW(110) layers were individually fully relaxed for 
i = 1 to 4 and j = to 3; j = refers to the uncovered 
iW(110). The conventional unit cell of the 3U/7W(110) 
system is highlighted, as shown in Fig. [1] (c). The prim- 
itive unit cell of the highlighted portion of the relaxed 
slab, which crystallizes in the Cmmm space group, con- 
stitutes our superlattice surface structure. 



IV. RELAXED INTERLAYER SPACINGS 

Vertical IL lattice spacings were calculated for the 
covered and uncovered W(110) surfaces and subsur- 
faces. The results after relaxations together with the 
experimental^ and theoretical^— results are presented 
in Tab. HI The relative distances in percentage were also 
given in this table. The relative changes for the two last 
shallow uranium layers were obtained by comparing with 
the bulk U-U vertical IL distance (ILD) of d^ (001) = 
2.493 A. The results are presented in columns 2 and 3 of 
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FIG. 2: (Color online) DOSs for (a) a-U compound, 
(b) bcc-W compound, (c) 1U/7W(110), (d) 1U/7W(110), 
(e) 2U/7W(110), (f) 3U/7W(110), (g) 1U/7W(110), (h) 
2U/7W(110), (i) 3U/7W(110). 



Tab.Q] For the first deposited uranium layer and all of the 
tungsten layers, as presented in columns 4 to 7 of Tab. U 
the relative changes were obtained by comparing with the 
bulk W-W vertical ILD of d^V(llO) = 2.246 A. For 
7 ML W(110) clean substrate, our result shows that the 
surface shallow vertical ILS d|^^_ W / 2 (110) is contracted 
to 2.140 Aby 4.7%, which is in agreement with the exper- 
imental valued of (-3.1±0.6)% and theoretical valued of 
-3.3%. Our calculated contraction d^_ lV2 (110) result 
for 5 ML W(110) clean system, 5.3%, is in agreement 
with the result of the others^, 4.1%, for the same sys- 
tem. Our 5W(110) and 7W(110) results show that the 
shallow ILS by adding one more bilayers is expanded from 
-5.3% for 5W(110) to -4.7% for 7W(110). The same trend 
was reported earlier^ as shown in Tab. HI However, this 
trend cannot be generalized to every layers, as a small 
expansion, which is not presented here, is occurred by go- 
ing from 7W(110) to 9W(110) in agreement with Ref^. 
The accuracy of the experimental results, as represented 
in Tab. U is about 1.0% for deeper layers of tungsten 
which are given in the last two columns of Tab. HJ For 
the later layers our results in agreement slightly with the 
theoretical results of the others^—, predict small re- 
laxations near the range of the experimental accuracy. 
The small relative changes can be then nearly neglected 
for these deeper layers of W(110). Now, let us grad- 
ually add the three uranium adatoms on the 7W(110) 
substrate layer by layer. By adding the first ML of ura- 
nium on the 7W(110) substrate (ultrathin 1U/7W(110) 
film), the result, as presented in Tab. HI shows that the 
vertical ILS at the interface, dfj{_ wi , is considerably ex- 
panded by 14% to 2.562 A when compared with the ver- 
tical ILS of the bec-bulk along the [110] direction, i.e., 
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TABLE I: Vertical ILDs for fully relaxed covered and uncovered tungsten substrate surfaces and subsurfaces in A. The percent- 
ages of the relative differences given in columns 2 and 3 of the table were obtained by comparing with the bulk U-U vertical 
ILD of d^"i fc [/ (001) = 2.493 A, while for the remained columns 4, 5, 6, and 7, they were obtained by comparing with the bulk 
W-W vertical ILD of d&r!V(H0) = 2.246 A. 



System 


jsur lsur 
a U3-U2 a U2-Ul 


jsur 

a Ul-Wl 


jsur 

a Wl-W2 


j sur 

a W2-W3 


jsur 

a W3~W4 


1U/7W(110) 




2.562 (14.0%) 


2.202 (-1.9%) 


2.196 (-2.2%) 


2.208 (-1.7%) 


2U/7W(110) 


2.018 (-19.1%) 


2.397 (6.7%) 


2.116 (-5.8%) 


2.114 (-5.9%) 


2.114 (-5.9%) 


3U/7W(110) 


2.102 (-15.7%) 2.164 (-13.2%) 


2.432 (8.3%) 


2.149 (-4.3%) 


2.143 (-4.6%) 


2.140 (-4.7%) 


5W(110) 






2.124 (-5.3%) 


2.226 (-0.9%) 




a 5W(110) 






(-4.1%) 


(-0.2%) 




b 5W(110) 






2.173 (-4.1%) 


2.258 (-0.4%) 




7W(110) 






2.140 (-4.7%) 


2.208 (-1.7%) 


2.219 (-1.2%) 


a 7W(110) 






(-3.3%) 


(-0.1%) 


(-0.4%) 


C 12W(110) 






(-4.7%) 


(0.2%) 


(-0.6%) 


n W(110) Ba=p ' 






(-3.1%±0.6%) 


(0.0%±0.9%) 


(0.0%±1.0%) 



"Reference^; The theoretical data were calculated within the APW method as embodied in the WIEN93 code using the 
PZ-LDA81. The experimental data were obtained from LEED analysis. 

6 Reference2£; The data were calculated within the FP-APW method as embodied in the WIEN97 code using the GGA92. 
c Reference-i; The data were calculated within the pseudopotential method as embodied in the Vasp code using the GGA92. 



d^^H/(H0) = 2.246 A. Because of the large expansion, 
one may expect to observe significant changes at the in- 
terface in the electronic structures and result in 
the sensitive physical properties, as shall be discussed 
in Sec. |V] and subsections therein . Conversely, for the 
2U/7W(110) system a dramatic contraction, 19.1%, is 
observed in Tab. Ufor the outer vertical ILS of d^ 2 _ ui 
= 2.018 A compared to the bulk of uranium ILS along 
the [001] direction, i.e., d^[ fc [/ (001) = 2.493 A. The con- 
traction gives rise to an enhancement of the direct 5f-5f 
interactions between two adjacent uranium atoms along 
the z-Cartesian direction. Consequently, one may ex- 
pect to observe more broadened density of states (DOS) 
for the 2U/7W(110) system compared to the 1U/W(110) 
system, as shall be shown in Figs. [5] By adding the 
third uranium layer to the last 2U/7W(110) slab and 
formation of the 3U/7W(110) system, the new outer 
ILD, df^_ U2 = 2.102 A, is slightly expanded compared 
to the last outer ILS of d s ^_ m = 2.018 A. But, the 
later distance of djj 3 r _ U2 is still contracted by 15.7% 
compared to the d^(001) = 2.493 A and by 19.8% 
compared to the outer ILS of dfj l [_ wi = 2.562 A for 
the ultrathin 1U/7W(110) film. Therefore, one may ex- 
pect to observe similar behavior for both of 3U/7W(110) 
and 2U/7W(110) systems. But these 3U/7W(110) and 
2U/7W(110) systems may be expected to behave differ- 
ently compared to the ultrathin 1U/7W(110) film, as 
shall be shown in Sees. IV Bl IV Dl and IV El As can be 
seen from columns 4 to 7 of Tab. [J our result shows that 
the ILSs oscillate versus uranium film thickness. For in- 
stance, the d^"[_ W2 first decreases from 1U/7W(110) to 
2U/7W(110), while it then increases from 2U/7W(110) 



to 3U/7W(110). The later oscillation can be also seen 
from the dff{_ wi , d$ r 2 _ W3 , and d$ r 3 _ W4 columns as 
well as the aforesaid dK_wo. 



V. MAGNETIC PROPERTIES AND 
ELECTRONIC STRUCTURES 



A. Bulks of U and W 



We calculated the MMs and the DOS for the or- 
thorhombic a-U. The DOS showed in Fig. [5] (a) indicates 
that the system is nonmagnetic in its bulk form. This 
result is in excellent agreement with the previous results 
of the other s 7 ' 11 within different methods. It was shown 
that the bulk of a-U compound could be considered as 
a nonmagnetic system using FP-APW method in Refj^i 
and using FP-LMTO method in Refi 7 -. The later group 7 , 
calculated spin and orbital MMs per uranium atom to be 
0.1 /xb and -0.22 fiB for the hypothetical hcp-U in its bulk 
form, respectively. This shows that the structure of the 
uranium can affect the magnetic properties of the ura- 
nium based systems. Thus, it may be important to take 
an actual structure into consideration for the uranium 
as the overlayer deposited on the tungsten substrate. We 
also made sure about the magnetism of the tungsten sub- 
strate in its bulk form. The bulk of tungsten is also found 
to be nonmagnetic, as may be seen from Fig.[2](b). 
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FIG. 3: (Color online) Freestanding 7 ML q-U(OOI). 



TABLE II: Spin MM inside the Muffin-tin sphere (MT), spin 
MM in the interstitial region (Int), total spin MM per U 
atom(SP), orbital MM (Orb) and total MM (Tot) in fj, B 
within GGA+SO for the freestanding a-U(OOl) TFs from 1 
ML to 7 ML. For more than 1 layer the total spin MMs per 
atom were calculated as MT+~Int. 





1 layer 


3 layers 


5 layers 


7 layers 


MT 


2.132 


0.867 


0.864 


0.865 


Int. 


0.754 


0.412 


0.409 


0.411 


SP 


2.886 


1.073 


1.068 


1.070 


Orb 


-0.812 


-0.372 


-0.372 


-0.373 


Tot 


2.074 


0.701 


0.696 


0.697 



B. Freestanding U Surfaces 

To reinforce our approach and validate its predictions 
on the actual a-U/W(110) case, it may be expedient to 
first try to reproduce the freestanding results. The free- 
standing results can be directly compared with the theo- 
retical results of the others. We constructed a slab con- 
taining 7 ML U, as illustrated in Fig. [31 to simulate the 



freestanding a-U. We increased the vacuum to reach the 
optimized value of 20 bohr, i.e. ^5.29x2 = 10.58A, as 
shown in Fig. |3l We made sure that total energy, sur- 
face energy and work function of the system were not 
changed by increasing the amount of vacuum. We then 
fully relaxed the uranium layers to minimize the exerted 
forces on each atoms, as mentioned in Sec. HH We cal- 
culated spin MM inside the Muffin-tin sphere, total spin 
MM, orbital MM inside the Muffin-tin sphere, interstitial 
and total MMs within the GGA+SO for the freestanding 
a-U surfaces. Our results are presented in Tab. HU The 
total spin MM can be obtained by adding the spin MM 
inside the Muffin-tin sphere to its interstitial contribu- 
tion. Stojic and collaborators^ calculated spin MM in- 
side the Muffin-tin sphere, orbital MM inside the Muffin- 
tin sphere, interstitial, and total MMs per U atom to be 
2.11 (0.84), -0.83 (-0.38), 0.74 (0.39), and 2.02 (0.66) in 
Hb for 1 ML (3 ML) of the hypothetical freestanding a-U 
TFs using the FP-APW method within the GGA+SO, 
respectively. Their results^ were well converged and did 
not show significant changes by adding more uranium 
ovcrlaycrs. A comparison shows that our freestanding re- 
sults, as presented in Tab.lLTl are in good agreement with 
the results given in Ref The small differences between 
our results and the results reported in Ref ^ could be due 
to the fact that they-ii did not relax their simulated free- 
standing a-U TFs and utilized the unrelaxed experimen- 
tal geometry. The agreement authenticates the validity 
of our freestanding results. Spin and orbital MMs per 
uranium atom were calculated in ReLi to be ~0.2 hb 
and ^-0.44 fis for the hypothetical freestanding hcp-U 
surface, respectively. Our result shows that the spin and 
orbital MMs of the first uranium over layer are aligned an- 
tiparallel in accordance with the above results reported 
in Refsiii. Spin MM and the magnitude of the orbital 
MM of a-U freestanding TFs are larger than those of the 
hcp-U freestanding TFs. This is in the case that, as dis- 
cussed in Sec. IV Al the hcp-U shows magnetic ordering in 
its bulk form, whereas the bulk of a-U is not a magnetic 
system. Hence, the structure of uranium can cause to 
reverse the magnetic properties from its bulk to its sur- 
faces and inversely from its surfaces to its bulk. Similar 
results were recently reported^ for the isostructural a- 
and 7-Ce systems. 



C. q-U/W(110) Thin Films 

Spin MMs per atom were calculated for the U and W 
in the jU/7W(110) system with j = 1, 2, 3 using the LDA 
and GGA excluding and including SO coupling. Orbital 
polarizations were included only for the U atoms in a 
non-self-consistent manner. The non-self-consistent pro- 
cedure seems to be adequate for predicting orbital mo- 
ment, because 5f U electrons in this system behave as 
band-like electrons. To obtain sclf-consistently the or- 
bital MMs of the U overlayers, one may include suit- 
able correlations among 5f-electrons, e.g., by performing 
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TABLE III: Spin MM inside the Muffin-tin sphere of 
Uj atom (MT(Uj) for j=l to j ma x), and of Wi atom 
(MT(Wi) for i=l to 4), spin MM in the interstitial re- 
gion (Int.), total spin MM of the primitive unit cell 
(SP), orbital MM inside the Muffin-tin sphere of Uj atom 
(ORB(Uj) for j=l to jmax), and total MM per primi- 
tive unit cell (TOT) in fi B in covered jU/7W(110) for j 
= 1 to jmax within the PW-LDA and PBE-GGA with 
and without SO coupling. The total spin moment is 
SP=2x(£j™r MT ( U i)+Et=i MT (Wi))+Int. Total MM per 

primitive unit cell (TOT) is T0T=SP+2x Ej=r 0RB ( U j)- 
The jmax index varies from 1 in 1U/7W to 3 in 3/7W(110). 
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LDA+U calculations. Here, instead, following the exper- 
imental results^— we assumed that the 5f electrons were 
not localized in U/W(110) TFs, as notified in SecHH The 
results together with the total MM per primitive unit cell 
are presented in Tab. lIIII One may observe some palpable 
trends which can be immediately deduced from Tab. IIIII 
The magnitude of our calculated MMs within the GGA 
is larger than that of the LDA. We can come to the same 
conclusion by comparing the result of the GGA+SO with 
that of the LDA+SO. This implies that the SO coupling 
does not affect the trend. However, SO gives rise to re- 
duce the MMs as can be observed by comparing the result 
of LDA (GGA) with that of the LDA+SO (GGA+SO), 
as presented in Tab. IIIII The inclusion of the SO cou- 
pling makes the up-DOS more symmetric with respect 
to the down-DOS compared to the SO excluded case 
(see Figs. [2] (d), and (g)). Thereby, the difference be- 
tween up- and down-DOSs becomes smaller by including 
SO coupling due to the more symmetric up- and down- 
DOSs. The more symmetric up- and down-DOSs may 
result in less MM. Therefore, one may expect to obtain 
smaller MM by including SO coupling for this system. 
The reduction of MM from GGA to GGA+SO is com- 
parable with the reduction of MM from GGA to LDA, 
as shown in Tab. [TTTJ Therefore, the results of GGA+SO 
and LDA seem to be approximately closer to each other 
than the results of the other exchange-correlation func- 
tionals (see Tab. IIII[) . Our surface MM calculation does 
not reveal any considerable MMs for Wi with i = 1 to 
4. Total up-DOS of WI cancels total down-DOS of WI 
(sec Fig. [21(c)). Therefore, our result, in agreement with 
Refj2£, shows that the tungsten even at its surfaces is 
not capable of producing magnetism. Our result pre- 
dicts spin MM at the first layer of uranium (Ul) for the 
ultrathin 1U/7W(110) system, as shown in Tab.|HU Or- 
bital MM is also predicted by our calculations including 
SO coupling in 1U/7W(110). The total- and f-DOSs for 
Ul in ultrathin 1U/7W(110) film, as shown in Fig. [1(d), 
corroborate that the nonmagnetic bulk a-U becomes a 
magnetic layer when its first ultrathin layer is deposited 
on the nonmagnetic W(110) substrate. The source of 
the imposed MM in the first layer of 1U/7W(110) may 
be attributed to the vertical ILD at the interface, as 
presented in Tab. [J] The result shows that the Ul-Wl 
distance, d^J l [_ wl = 2.256 A, is expanded by 14% com- 
pared to d^V(lOO) = 2.246 A, and by 2.8% compared 
to d^i fc y(100) = 2.493 A. More importantly, there is 
no any other uranium layer over the first deposited U 
layer in the case of 1U/7W(110). Thus, one does not 
account for the strong direct 5f-5f interactions perpen- 
dicular to the surface. This gives rise to a more localized 
DOS near the Ep, as shown in Fig. [2] (d). Thereby, 
there is no surprise in observing magnetism in the ul- 
trathin 1U/7W or in the freestanding 1ML U system. 
Our result, as shown in Tab. IIIII shows that by adding 
the second U layer to the 1U/7W(110) or by adding the 
third layer to the 2U/7W(110), all the components of 
the total MMs are substantially reduced. These sensi- 
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TABLE IV: Vertical ILDs for fully relaxed 3 ML, 5 ML and 
7 ML freestanding a-U in A. 



System 


j sur 
a U2-Ul 


isur jsur isur 
a U3-U2 a U4-U3 a U'2-Ul 


a V'3-U'2 a U'4,-U' 


3U(001) 


2.404 


2.405 




5U(001) 


2.340 


2.275 2.342 


2.278 


7U(001) 


2.325 


2.261 2.187 2.326 


2.263 2.188 



ble reductions in the MM may be attributed to the ti- 
ll ILSs in the 2U/7W(110) or 3U/7W(110), which are 
smaller than the Hill limi^. According to the Hill crite- 
ria, uranium based compounds with the U-U interatomic 
spacings smaller than 3.4-3.6 A are typically nonmag- 
netic with itinerant 5f states^. Therefore, 2U/7W(110) 
with df^_ m = 2.018 A and 3U/7W(110) with df%_ m = 
2.164 A and df%_ U2 = 2.102 A, would be almost non- 
magnetic based on the Hill limit. Our LDA+SO and 
nearly LDA as well as GGA+SO results show that the 
2U/W(110) and 3U/W(110) systems are not magnetic, as 
shown in Tab. [Hi in agreement with the Hill prediction. 
Our calculated DOSs, as shown in the Figs. [2] (e), (f), 
(h) and (i), authenticate that the 2U/7W(110) and/or 
3U/7W(110) TFs become nonmagnetic compared to the 
magnetic 1U/7W(110) ultrathin film, as shown in the 
Figs. [5] (d) and (g), by absorbing additional uranium 
over layers. 



D. Comparison of freestanding a-U(OOl) and 
a-U/W(110) Thin Films 

The calculated spin MM inside the Muffin-tin sphere 
of the 1U/7W(110), 1.402 (i B , as given in Tab. [THl is 
smaller than the value of 2.1 Ms for the freestanding 1 
ML a-U, as presented in Tab. Ullor in RefAi. If we com- 
pare the interstitial MM of the 1U/7W(110), 1.088 Ms, 
as given in Tab. IIIIl with that of 0.75 or 0.74 Ms for 
the freestanding 1 ML a-U, as presented in Tab. [XT] or 
Rebii, we should first divide the one for 1U/7W(110) by 
2, i.e., 1.088/2 = 0.544 Ms- In our 1U/7W case for 1 
ML U there are two U atoms, which is not the case for 
the 1 ML of the freestanding a-U system. For the later 
case, 1 ML freestanding a-U, there is only one U atom 
in the supercell. The comparison shows that the value of 
0.74 hb for the freestanding a-U is larger than the cal- 
culated value of 0.544 mb for the 1 ML a-U deposited on 
the W(110) substrate. The difference may be due to the 
proximity effect of the substrate. All of the atoms in the 
superlattice structure contribute to the interstitial MM, 
and there at the interstitial region their contributions are 
mixed with each other. However, similar to comparing 
total MM it is enough to divide the interstitial moments 
by 2, as the W atoms do not contribute substantially in 
the total MM. The comparison shows that the GGA+SO 
total MM per U atom, -1.402+1.088/2-0.446=1.5 Ms, as 



given in Tab. IIIIl regardless of the ignorablc W atoms' 
contributions, is less than the freestanding 1 ML total 
MM of 2.074 or 2.02 mb, as presented in Tab. HI or in 
Refiii. This may demonstrate the effect of substrate. To 
reconfirm this fact, we can compare the orbital MMs as 
well. Even though the signs of the orbital MM for the 
1U/7W(110) and the freestanding a-U are the same, the 
magnitude of the GGA+SO result for the 1U/7W(110), 
0.446 /is, is half of the value for 1 ML freestanding a-U, 
0.8 Ms- This shows that all the components of the to- 
tal MM for the 1U/7W(110) are reduced when compared 
with those of the freestanding 1 ML a-U. The reduction 
may come from the nearness effect of substrate on the 
TF magnetism. In addition to the discussed compari- 
son between 1 ML freestanding U and the 1U/7W(110), 
one may compare the trend of the MMs for the free- 
standing a-U by adding more uranium layers with that 
of the aU/7W(110). A significant reduction in total MM 
from -2.0 mb for 1 ML to -0.7 /i B for 3 ML can be 
observed in the freestanding a-U TFs, as presented in 
Tab. Ullor in Ref^i. This shows that the number of layers 
can be of significant importance. However, we still ob- 
serve the nonzero total MM of —0.7 mb for 3 and larger 
ML of the freestanding a-U, as presented in Tab. |H] or 
in Ref^i. Based on the hypothetical freestanding a-U 
model in agreement with RefAi, one predicts the surface 
magnetism. Our GGA+SO result predicts small MMs 
of 0.057 mb and -0.047 Ms P er primitive unit cell for 
the 2U/7W(110) containing 4 U atoms and 3U/7W(110) 
containing 6 U atoms, respectively. Therefore, the maxi- 
mum total MM per U atom that can be predicted by our 
GGA+SO result is only 0.0014 mb for 2U/7W(110) and 
-0.0078 /is for 3U/7W(110) TFs. The later tiny value of 
0.0014 mb per U for 3U/7W(110) is much smaller than 
the value of — 0.7ms per U for the freestanding 3 ML a-U. 
In essence, our result based on the actual a-U/W(110) 
system predicts that the 3U/W(110) is not a magnetic 
system, while the result for the freestanding case shows 
that the 3 ML a-U is a magnetic system. In order to find 
the physics behind the later result, we calculated the ver- 
tical IL lattice spacings for the freestanding surfaces, as 
shown in Fig. [3J The results after relaxations were pre- 
sented in Tab. HVl The vertical ILD before relaxation, as 
shown in Fig.[T](b), is the distance between the uranium 
layers of the a-U in its bulk form along the [001] direc- 
tion, i.e., d^(001) = 2.493 A. Although the result, as 
can be seen from Tab. IIV1 shows that the vertical ILSs 
for the relaxed freestanding surfaces, as shown in Fig. [31 
are contracted by relaxation compared to the unrelaxed 
experimental geometry, as shown Fig. [1] (b), the vertical 
ILSs are still larger than those for the a-U/W(110) TFs, 
as can be seen from Tab. [U For instance, the vertical 
ILSs are 2.404 and 2.405 A for the freestanding 3 ML 
a-U(OOl) TFs, which are larger than those of 2.102 and 
2.164 A for the 3U/7W(110), as can be seen from Tabs.[fl 
and IIVI The smaller U-U spacings, the stronger U-U in- 
teractions among their 5f electrons. The stronger 5f U-5f 
U interaction gives rise to more broadened 5f DOS, which 



8 



TABLE V: Total MM in \i B per primitive unit cell within the 
PBE-GGA for three cases when the Ul is (i) in its relaxed 
equilibrium position, dy" r -wi =2.526 A, (ii) moved upwards 
by 0.50 A, d£f_ m =3.026 A, (iii) moved upwards by 1.00 A, 
df/i-wi =3.526 A, with respect to the Wl, while the other 
ILDs given in Tab. [T] were kept fixed. 





duT-wi+0 A 


d^T-wn+0.5 A 


dhJ-wx+1 A 


1U/7W(110) 


4.56 


5.39 


5.54 


2U/7W(110) 


0.45 


0.53 


0.57 


3U/7W(110) 


0.28 


0.33 


0.35 



results in less localization. Therefore, one may expect to 
obtain less MM for the 3U/7W(110) system compared to 
the freestanding 3 ML a-U(OOl) TFs. Our result verifies 
such an expectation, as shown in Tabs. Q] and IIVI Our 
result, as shown in Tab. IIVI shows that by adding more 
layers, the vertical IL U-U spacings are reduced for the 
freestanding U TFs. For example, the distance between 
Ul and U2 is reduced from 2.404 A for 3 ML freestand- 
ing U(001) to 2.340 A for 5 ML freestanding U(001), 
as shown in Tab. IIVI This reduction may explain why 
the MM is reduced from 0.701 fiB for 3 ML freestanding 
U(001) to 0.696 [i B for 5 ML freestanding U(001). 



E. Proximity and lattice mismatch effects of 
W(110) substrate 

In Sec. IV Dl we related the larger MM of the free- 
standing a-U(OOl) than that of the a-U/W(110) to the 
contracted vertical IL U-U distances of the a-U/W(110) 
compared to that of the freestanding a-U(OOl). How- 
ever, there we did not discuss the role of the tungsten 
substrate. The effect of tungsten substrate may become 
more important when we compare the MM of the free- 
standing 1 ML a-U(OOl) with that of 1U/7W(110). The 
MM per U atom within the GGA+SO for the freestand- 
ing 1 ML a-U(OOl), 2.074 fi B , is larger than that for the 
1U/7W(110), 2.988/2 = 1.494 fj, B , as shown in Tabs. HI 
and IIIII This is in the case that there is no other uranium 
layer over the first uranium layer for both 1 ML a-U(OOl) 
and 1U/7W(110). Therefore, for these one ML U cases 
there are not any vertical U-U interactions. Thus, the 
vertical U-U interactions could not be the case for the 
contraction and as a result for the reduction of the MM 
of the 1U/7W(110) compared to that of the 1 ML a- 
U(001). In this section, we show that, in addition to the 
vertical IL U-U spacings, the in-plane U-U spacings may 
also affect the MM. In order to achieve this goal, we first 
discuss the proximity effect to elucidate the role of tung- 
sten on its first uranium adatoms through the vertical 
W-U interaction at the interface of the 1U/7W(110) sys- 
tem. We then try to show the effect of tungsten through 
its imposed lattice mismatch on the in-plane U-U spac- 
ings. We discuss that not only the vertical IL U-U in- 



teractions, but also the in-plane U-U interactions can in- 
fluence the MM. To quantitatively elucidate the effect 
of the distance between the substrate and the film, we 
calculated, as presented in Tab. [Vj the total MM per 
primitive unit cell within the GGA when the Ul is ver- 
tically displaced 0.5 A and 1.0 A upwards with respect 
to its relaxed equilibrium position at df] l {_ wi = 2.526 A. 
In the displacement procedure, it is assumed that apart 
from the dfj{_ wl distance all of the other distances are 
kept fixed with no changes in their relaxed values as given 
in Tab.[U The result, as presented in Tab. IVl shows that 
the total MMs increase as the Ul together with its fixed 
overlayers, if any, is moved upwards from its equilibrium 
position with respect to the fixed nonmagnetic 7W(110) 
substrate. The total MMs increase by ^15% when the 
layer is pulled up by 0.5 A. The percentage docs not 
change significantly (more or less ^21%) when the layer 
is pulled up by 1.0 A compared to the equilibrium posi- 
tion. Total MM was found±i to be 3.27 [i B per U for 1 ML 
a-U within GGA. Our GGA result for the 1U/7W(110) 
is 4.56 \ib per primitive unit cell and thereby 4.56/2 = 
2.28 fiB per U atom when Ul is placed at its equilib- 
rium position. It increases from 2.28 fiB to 5.54/2 = 
2.77 [lb as the overlayer is pulled away by 1 A from 
the tungsten substrate. Thus the proximity effect of the 
substrate can cause our GGA result of 2.77 fiB for the 
1U/7W(110) to be closer to the GGA result 1 ! of 3.27 
\xb for the freestanding system compared to the value 
of 2.28 fiB for 1U/7W(110) before its uranium overlayer 
is moved up. However, our result is still 15% smaller 
than the freestanding result. The source of the later dif- 
ference may not be only due to the proximity effect of 
the substrate, as further increase in vertical spacing of 
dui-wi d° es n °t affect the result substantially. Indeed, 
it may originate from the in-plane lattice mismatch ef- 
fect of the substrate. The in-plane lattice mismatches 
can cause the lattice parameters of a-U, as given in the 
text and Figs.[TJ to expand by 36.7% along a axis and to 
contract by 45.1% along b axis. Therefore, the difference 
may not be considered so large by accounting the large 
in-plane lattice mismatches. More precisely, the in-plane 
U-U distance between n earest neighbors for the case of 
the freestanding a-U is V / (2M2) 2 + (5.780) 2 /2 ~ 3.22 A, 
which is larger than that for the a-U/7W(110) system, 
i.e., ^(4.492)2 + (3.176) 2 /2 ~ 2.75 A. For deriving the 
former value, we have used the lattice parameters of the 
a-U, as shown in Fig. [T] (b), while for deriving the later 
value we have used the lattice parameters of the W(110) 
substrate, as shown in Fig. [I] (a). Although both of these 
distances are smaller than the Hill limit, the in-plane U- 
U distance for the freestanding is larger and as a result 
closer to the Hill limit than that for the a-U/W(110) 
system. Consequently, it seems to be reasonable to find 
larger MM for the freestanding case compared to the MM 
of the supported case by the tungsten substrate. In ad- 
dition to the in-plane lattice mismatches and the prox- 
imity effects of the substrate, vertical lattice mismatch 
of the substrate can be also of significant importance in 
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FIG. 4: (Color online) Band structures of 1U/7W(110), 2U/7W(110) and 3U/7W(110) systems. The light lines indicate the 
total and the dark lines indicate the uranium f-orbital bands. 



this 5f system. In contrast to the ultrathin 1U/7W(110) 
film, the 2U/7W(1I0) and 3U/W(II0) TFs, as shown in 
Tab. IIV1 do not manifest considerable magnetism even 
by moving upwards the Ul overlayer from the nonmag- 
netic substrate. This once more shows that the proxim- 
ity effect of the substrate cannot be the only case here, 
since it has been already showr>ii that the freestanding 
3 ML a-U is a magnetic system. As discussed in Sec. lIVl 
the distance between uranium overlayers in 2U/7W(110) 
system is remarkably contracted by 19.1%, as shown in 
Tab. HI compared to the bulk of a-U. The contraction 
for the 3U/7W(110) is 15.7% between U3 and U2 and 
13.2% between U2 and Ul. The contractions can cause a 
stronger direct interaction between 5f states of two ver- 
tically adjacent uranium atoms. The stronger 5f-5f in- 
teraction causes the valence DOS to be more broadened. 
Total- or 5f-DOS of Ul, as shown in Figs. [2] (d) and (g), 
for U1/7W(110) shows a nearly narrow sharp peak in 
the vicinity of Ep. However, by adding more uranium 
overlayers the sharp narrow peak is more broadened, as 
can be seen from Figs. [2] (e), (f), (h), (i) compared to 
Fig. [2] (d) and (g). Hence, one may anticipate to ob- 
serve strong reduction in magnetism from 1U/7W(110) 
to 2U/7W(110) or 3U/7W(110) in complete accordance 
with our result presented in Tabs. IIIII and [V] The com- 
parison between Fig. [2] (c) or (f) with Fig. [2] (d) or be- 
tween Fig. [5] (h) or Fig. [2] (i) with Fig. [5] (g) corroborate 
the anticipation, since due to the direct interactions the 
DOSs are more symmetric and broadened. If we compare 
total- and f-DOSs of Ul with each other for the ultrathin 
Ul/W(110) film, as shown in Fig. [2] (d), we see that, 
in complete accord with experiments^^, most of the 
total-DOS is dominated by the 5f states in the vicinity 
of the Ep. In order to authenticate the above discussion 
we have calculated the band structures of 1U/7W(110), 
2U/7W(110) and 3U/7W(110) systems. Our results are 
shown in Fig. [4] The calculated band structures show 
that the 5f U band character dominates in the vicinity 
of Ep in agreement with experiments^ and our calcu- 
lated DOS. The comparison of band structures between 



1U/W(110) and 2U/W(110) TFs shows that the degener- 
acy of 5f band is lifted by the vertical U-U interactions at 
-1 cV. The removed degeneracy causes the band struc- 
tures of the 2U/7W(110) to be more dominated in the 
vicinity of the Ep by the 5f character of the uranium 
compared to that of the ultraTF of 1U/W(110). The 
number of lifted 5f bands increases in 3U/W(110), as 
shown in Fig. 21 as the third uranium layer is deposited 
on the 2U/W(110) TFs. Such an increase in the num- 
ber of splitted 5f bands in the 3U/W(110) can be taken 
as more interactions among 5f electrons of the uranium 
atoms in this system. This shows that strong interactions 
among 5f U electrons in uranium TFs in the presence of 
tungsten substrate may not allow the 3U/W(110) system 
to manifest magnetization at its surfaces. In conclusion 
the tungsten substrate can cause to destroy the appear- 
ance of magnetization on the hypothetical freestanding 
a-U(OOl) TFs. 



VI. SUMMARY AND OUTLOOK 

In this study, wc considered the a-structure for the ura- 
nium overlayers based on the recent experimental obser- 
vations. We transformed the surface of the orthorhombic 
a- uranium along the [001] direction, and found it to be 
similar to that of the bec-tungsten substrate along the 
[110] direction. The surface structures of the W(110) 
and the transformed a-U(OOl) are both fcc-like, as dis- 
cussed in the text. We notice that the a-U might be 
an ideal adsorbate for the W(110) tungsten taking ver- 
tical and in-plane lattice mismatches into consideration. 
This is in agreement with the recent experiment. Wc 
then simulated the a-U/W(110) TFs and calculated its 
surface and electronic properties. The calculations were 
performed within the DFT employing the accurate FP- 
APW+lo method as embodied in the reliable WIEN2k 
code. Our result shows the proximity effect of the non- 
magnetic tungsten substrate on the magnetism of ura- 
nium overlayers. In addition to the proximity effect, 
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the imposed vertical and in-plane lattice mismatches of 
the W(110) substrate were found to be of significant 
importance. Our result predicts surface magnetism in 
the ultrathin 1U/7W(110) film. This result is in agree- 
ment with that of the freestanding 1ML a-U(OOl). By 
adding more uranium layers on the ultrathin film, we 
found that the imposed surface magnetism disappears. 
This is in the case that for the freestanding a-U(OOl) 
TFs, the magnetism is preserved but its magnitude de- 
creases by adding one more uranium layer. The disap- 
pearance of the magnetism in 2U/7W(110) is attributed 
to the contraction of the uranium ILS, which can be 
originated from the W(110) substrate. The contraction 
can give rise to stronger direct 5f-5f interactions. Our 
calculation shows that the 5f U DOS is more broad- 
ened in 2U/7W(110) or 3U/7W(110) compared to the 
ultrathin 1U/7W(110) film. The broadening of the 5f 
DOS corroborates that the 5f states is more band-like 
for 3U/7W(110). This is in agreement with the exper- 
imentally founded dispersive behavior for the 5f states 



in U/W(110) TFs. Our result in agreement with experi- 
ment shows that the 5f character predominates the total 
DOS in U/W(110) TFs. According to our results, we 
anticipate further theoretical and experimental elabora- 
tions on the U/W(110) TFs taking the ev-structure for 
the adsorbate uranium into consideration. For example, 
various orientations of the a-U ovcrlayers may improve 
the results, provided that less lattice mismatches can be 
found with the tungsten substrate. Theoretical investiga- 
tion of hcp-U/W(110) TFs may also provide a significant 
contribution to the filed. 
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